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Abstract 


PGA (ProGram Algebra) is an algebra of programs which concerns 
programs in their simplest form: sequences of instructions. Molecular 
dynamics is a simple model of computation developed in the setting of 
PGA, which bears on the use of dynamic data structures in program- 
ming. We consider the programming of an interpreter for a program 
notation that is close to existing assembly languages using PGA with 
the primitives of molecular dynamics as basic instructions. It happens 
that, although primarily meant for explaining programming language 
features relating to the use of dynamic data structures, the collection 
of primitives of molecular dynamics in itself is suited to our program- 
ming wants. 


1 Introduction 


In this paper, we consider the programming of an interpreter for a program 
notation that is close to existing assembly languages using PGA (ProGram 
Algebra). With that we carry on the line of research with which a start was 
made in [3]. The object pursued with that line of research is the development 
of a theoretical understanding of possible forms of sequential programs, 
starting from the simplest form of sequential programs, and associated ways 
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of programming. The view is taken that sequential programs in the simplest 
form are sequences of instructions. PGA, an algebra of programs in which 
programs are looked upon as sequences of instructions, is taken for the 
basis of the development aimed at. Until now, the research was rather 
focussed on the development of a theoretical understanding of possible forms 
of sequential programs. The work presented in the current paper is primarily 
concerned with programming using the simplest form of sequential programs 
and molecular dynamics. 

For the programming of the interpreter, we use PGA with the primi- 
tives of molecular dynamics as basic instructions. Molecular dynamics is a 
simple model of computation bearing on the use of dynamic data structures 
in programming. In this model, states of computations resemble collections 
of molecules composed of atoms and computations take place by means 
of actions which transform the structure of molecules like in chemical re- 
actions. Molecular dynamics has been developed mainly in the setting of 
PGA. The model introduced in the current paper is a small expansion of 
the model that was first described informally in [1] and formally in [6]. 

In the line of research carried on, the view is taken that the behaviours 
exhibited by sequential programs on execution are threads as considered 
in BTA (Basic Thread Algebra).? A thread proceeds by doing steps in a 
sequential fashion. A thread may do certain steps for having itself affected 
by some service or affecting some service. In the current paper, we will use a 
generalization of the use mechanism introduced in [8] and a complementary 
mechanism of that mechanism for such interactions between threads that 
are behaviours exhibited by sequential programs on execution and services 
that deal with steps that relate to molecular dynamics. A slightly different 
generalization of the use mechanism from [8] is introduced in [7] under the 
name action transforming thread-service composition. 

A hierarchy of program notations rooted in PGA is introduced in [3]. 
In the current paper, we consider the programming of an interpreter for one 
program notation that belongs to this hierarchy. The program notation in 
question, called PGLD (ProGramming Language D), is a very simple pro- 
gram notation which is close to existing assembly languages. The hierarchy 
also includes a program notation, called PGLS (ProGramming Language for 
Structured programming), that supports structured programming by offer- 


3In [3], basic thread algebra is introduced under the name basic polarized process 
algebra. Prompted by the development of thread algebra [8], which is a design on top of 
it, basic polarized process algebra has been renamed to basic thread algebra. 
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ing a rendering of conditional and loop constructs instead of (unstructured) 
jump instructions. Each PGLS program can be translated into a semanti- 
cally equivalent PGLD program. In [4], a variant of PGLD with indirect 
jump instructions is introduced. In the current paper, we show how to adapt 
the interpreter for PGLD to the presence of indirect jump instructions. 

This paper is organized as follows. First, we review BTA, PGA, and 
PGLD (Sections 2, 3, and 4). Next, we extend BTA with the mechanisms 
for interaction between threads and services used, introduce a state-based 
approach to describe services, and give a state-based description of services 
for molecular dynamics (Sections 5, 6, and 7). Following this, we give PGA 
programs for creating representations of PGLD programs by molecules and 
a PGA program for interpreting those representations (Sections 8 and 9). 
After that, we introduce the variant of PGLD with indirect jump instruc- 
tions and adapt the above-mentioned PGA programs to the presence of 
indirect jump instructions (Sections 10 and 11). Finally, we make some 
concluding remarks (Section 12). 


2 Basic Thread Algebra 


In this section, we review BTA, a form of process algebra which is tailored to 
the description of the behaviour of deterministic sequential programs under 
execution. The behaviours concerned are called threads. 

In BTA, it is assumed that there is a fixed but arbitrary finite set of 
basic actions A with tau ¢ A. We write Atay for AU {tau}. The members 
of Atay are referred to as actions. 

The intuition is that each basic action performed by a thread is taken 
as a command to be processed by a service provided by the execution envi- 
ronment of the thread. The processing of a command may involve a change 
of state of the service concerned. At completion of the processing of the 
command, the service produces a reply value. This reply is either T or F 
and is returned to the thread concerned. 

Although BTA is one-sorted, we make this sort explicit. The reason 
for this is that we will extend BTA with an additional sort in Section 5. 

The algebraic theory BTA has one sort: the sort T of threads. To build 
terms of sort T, BTA has the following constants and operators: 


e the deadlock constant D: T; 


e the termination constant S: T; 
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Table 1: Axiom of BTA 
xdtauby=adtaubse Tl 


e for each a € Atay, the binary postconditional composition operator 
_dab_:TxT-T. 


Terms of sort T are built as usual (see e.g. [13, 14]). Throughout the paper, 
we assume that there are infinitely many variables of sort T, including 
L,Y, 2. 

We use infix notation for postconditional composition. We introduce 
action prefizing as an abbreviation: ao p, where p is a term of sort T, 
abbreviates p Jab p. 

Let p and q be closed terms of sort T and a € Atay. Then pda gq will 
perform action a, and after that proceed as p if the processing of a leads 
to the reply T (called a positive reply), and proceed as q if the processing 
of a leads to the reply F (called a negative reply). The action tau plays a 
special role. It is a concrete internal action: performing tau will never lead 
to a state change and always lead to a positive reply, but notwithstanding 
all that its presence matters. 

BTA has only one axiom. This axiom is given in Table 1. Using 
the abbreviation introduced above, axiom Tl can be written as follows: 
x <itaub y=tauog. 

Each closed BTA term of sort T denotes a finite thread, i.e. a thread of 
which the length of the sequences of actions that it can perform is bounded. 
Guarded recursive specifications give rise to infinite threads. 

A guarded recursive specification over BTA is a set of recursion equa- 
tions E = {X =tx | X € V}, where V is a set of variables of sort T and 
each ty is a term of the form D, S ort dat’ with ¢ and t’ BTA terms of 
sort T that contain only variables from V. We write V(£) for the set of all 
variables that occur on the left-hand side of an equation in &. We are only 
interested in models of BTA in which guarded recursive specifications have 
unique solutions, such as the projective limit model of BTA presented in [2]. 
A thread that is the solution of a finite guarded recursive specification over 
BTA is called a finite-state thread. 

We extend BTA with guarded recursion by adding constants for solu- 
tions of guarded recursive specifications and axioms concerning these ad- 
ditional constants. For each guarded recursive specification E’) and each 
X € V(E), we add a constant of sort T standing for the unique solution 
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Table 2: Axioms for guarded recursion 
(X|E) =(tx|E) if X=tx ¢ E RDP 
E => X=(X|E) if X €V(E) RSP 


Table 3: Approximation induction principle 
An>o Tn(®) = Trly) > «= y AIP 


of E for X to the constants of BTA. The constant standing for the unique 
solution of FE for X is denoted by (X|E). Moreover, we add the axioms for 
guarded recursion given in Table 2 to BTA, where we write (tx|E) for tx 
with, for all Y € V(£), all occurrences of Y in tx replaced by (Y|E). In this 
table, X, tx and F stand for an arbitrary variable of sort T, an arbitrary 
BTA term of sort T and an arbitrary guarded recursive specification over 
BTA, respectively. Side conditions are added to restrict the variables, terms 
and guarded recursive specifications for which X, ty and EF stand. The 
equations (X|F) = (tx|E) for a fixed E express that the constants (X|E) 
make up a solution of E. The conditional equations E > X = (X|E) 
express that this solution is the only one. 

We will use the following abbreviation: a”, where a € Atay, abbreviates 
(X|{X =ao X}). 

We will write BTA+REC for BTA extended with the constants for 
solutions of guarded recursive specifications and axioms RDP and RSP. 

In [5], we show that the threads considered in BTA+REC can be viewed 
as processes that are definable over ACP [10]. 

Closed terms of sort T from the language of BTA+REC that denote 
the same infinite thread cannot always be proved equal by means of the 
axioms of BTA+REC. We introduce the approximation induction principle 
to remedy this. The approximation induction principle, AIP in short, is 
based on the view that two threads are identical if their approximations up 
to any finite depth are identical. The approximation up to depth n of a 
thread is obtained by cutting it off after performing a sequence of actions 
of length n. 

AIP is the infinitary conditional equation given in Table 3. Here, follow- 
ing [3], approximation of depth n is phrased in terms of a unary projection 
operator 7, :'T — T. The axioms for the projection operators are given in 
Table 4. In this table, a stands for an arbitrary member of Atay. 
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Table 4: Axioms for projection operators 


to(a) =D PO 
Tn4i(S) =S Pl 
Tn41(D) = D P2 
Tnii(@ Jab y) = m,(x) dal am,(y) P3 


We will write BTA+REC+AIP for BTA+REC extended with the pro- 
jection operators and the axioms from Tables 3 and 4. 


3 Program Algebra 


In this section, we review PGA, an algebra of sequential programs based on 
the idea that sequential programs are in essence sequences of instructions. 
PGA provides a program notation for finite-state threads. 

In PGA, it is assumed that there is a fixed but arbitrary finite set 2 of 
basic instructions. PGA has the following primitive instructions: 


e for each a € 2, a plain basic instruction a; 


for each a € 2, a positive test instruction +a; 


for each a € 2, a negative test instruction —a; 


for each 1 € N, a forward jump instruction #1; 
e a termination instruction !. 


We write 3 for the set of all primitive instructions. 

The intuition is that the execution of a basic instruction a may modify 
a state and produces T or F at its completion. In the case of a positive 
test instruction +a, basic instruction a is executed and execution proceeds 
with the next primitive instruction if T is produced and otherwise the next 
primitive instruction is skipped and execution proceeds with the primitive 
instruction following the skipped one. In the case where T is produced and 
there is not at least one subsequent primitive instruction and in the case 
where F is produced and there are not at least two subsequent primitive 
instructions, deadlock occurs. In the case of a negative test instruction 
—a, the role of the value produced is reversed. In the case of a plain basic 
instruction a, the value produced is disregarded: execution always proceeds 
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Table 5: Axioms of PGA 
(x;y);2=a;(y;z) PGAI 
(2) 2. = pe PGA2 
sy= a” PGA3 
zy)” =z; (y;x)” PGA4 


as if T is produced. The effect of a forward jump instruction #1 is that 
execution proceeds with the /-th next instruction of the program concerned. 
If 7 equals 0 or the /-th next instruction does not exist, then #/ results 
in deadlock. The effect of the termination instruction ! is that execution 
terminates. 

PGA has the following constants and operators: 


e for each u € J, an instruction constant wu; 
e the binary concatenation operator _;_; 
e the unary repetition operator _“. 


Terms are built as usual. Throughout the paper, we assume that there are 
infinitely many variables, including z, y, z. 

We use infix notation for concatenation and postfix notation for repe- 
tition. 

Closed PGA terms are considered to denote programs. The intuition is 
that a program is in essence a non-empty, finite or periodic infinite sequence 
of primitive instructions.4 These sequences are called single pass instruction 
sequences because PGA has been designed to enable single pass execution 
of instruction sequences: each instruction can be dropped after it has been 
executed. Programs are considered to be equal if they represent the same 
single pass instruction sequence. The axioms for instruction sequence equiv- 
alence are given in Table 5. In this table, n stands for an arbitrary natural 
number greater than 0. For each n > 0, the term x” is defined by induc- 
tion on n as follows: 2! = 2 and 2"t! = 2; 2". The unfolding equation 
xz’ = 2;2” is derivable. Each closed PGA term is derivably equal to a term 
in canonical form, i.e. a term of the form P or P; Q”, where P and Q are 
closed PGA terms that do not contain the repetition operator. 


4A periodic infinite sequence is an infinite sequence with only finitely many subse- 
quences. 
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Table 6: Defining equations for thread extraction operator 


al|=aoD #1| = D 
a;2|=aoc|z| #0;2|=D 
+a|=aoD #1; 2| = |2| 
t+a;2|=|z|dab|#2;2| |#l+2;uj=D 
—a|=aoD #1+2;u;2| = |#l+1;2| 
—a;2|=|#2;2|/dab lz] |!|=s 

fs o5 


Table 7: Rule for cyclic jump chains 
c= #0;y = |z|=D 


Table 8: Defining formulas for structural congruence predicate 
#n+13u.;...3;Un;#0 = #03 u1;...3;Un; #0 


#n+1sjur3...3Un;#m=H#mM+n4+13u1;...;tnj3#m 
(Fn +0+15ur5...5Un)® = (Al; ur; ..-5 un)” 
#m+nt+l+23ur3...3Unj3(13---5Um41)” & 

Pet LE ty FSi tig COT aS tra 
Ca 


Ly = yi AL2 = Yy2 > 115%. =M js y2Au1* Sy” 


Each closed PGA term is considered to denote a program of which the 
behaviour is a finite-state thread, taking the set 2 of basic instructions for 
the set A of actions. The thread extraction operator |_| assigns a thread to 
each program. The thread extraction operator is defined by the equations 
given in Table 6 (for a € &, 1 € N and u € 3) and the rule given in Table 7. 
This rule is expressed in terms of the structural congruence predicate _ = _, 
which is defined by the formulas given in Table 8 (for n,m,l € N and 
Wises lis Uiyont ig Opes) 

The equations given in Table 6 do not cover the case where there is a 
cyclic chain of forward jumps. Programs are structural congruent if they 
are the same after removing all chains of forward jumps in favour of single 
jumps. Because a cyclic chain of forward jumps corresponds to #0, the rule 
from Table 7 can be read as follows: if x starts with a cyclic chain of forward 
jumps, then || equals D. It is easy to see that the thread extraction operator 
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assigns the same thread to structurally congruent programs. Therefore, the 
rule from Table 7 can be replaced by the following generalization: « = y > 
|z| = lyl. 

Let E be a finite guarded recursive specification over BTA, and let Py 
be aclosed PGA term for each X € V(E). Let E’ be the set of equations that 
results from replacing in E all occurrences of X by |Px| for each X € V(E). 
If E’ can be obtained by applications of axioms PGA1—PGA4, the defining 
equations for the thread extraction operator and the rule for cyclic jump 
chains, then |Px| is the solution of E for X. Such a finite guarded recursive 
specification can always be found. Thus, the behaviour of each closed PGA 
term, is a thread that is definable by a finite guarded recursive specification 
over BTA. Moreover, each finite guarded recursive specification over BTA 
can be translated to a closed PGA term of which the behaviour is the 
solution of the finite guarded recursive specification concerned (cf. Section 4 
of [12]). 

Closed PGA terms are loosely called PGA programs. PGA programs in 
which the repetition operator does not occur are called finite PGA programs. 


4 The Program Notation PGLD 


In this section, we review a program notation which is rooted in PGA. 
This program notation, called PGLD, belongs to a hierarchy of program 
notations introduced in [3]. PGLD is close to existing assembly languages. 
It has absolute jump instructions and no explicit termination instruction. 

In PGLD, like in PGA, it is assumed that there is a fixed but arbitrary 
finite set of basic instructions 2. Again, the intuition is that the execution 
of a basic instruction a may modify a state and produces T or F at its 
completion. 

PGLD has the following primitive instructions: 


e for each a € 2, a plain basic instruction a; 
e for each a € 2, a positive test instruction +a; 
e for each a € A, a negative test instruction —a; 
e for each 1 € N, a direct absolute jump instruction ##1. 
PGLD programs have the form u1;...; uz, where ui,..., uz are primitive 


instructions of PGLD. 
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The plain basic instructions, the positive test instructions, and the 
negative test instructions are as in PGA. The effect of a direct absolute 
jump instruction ##I is that execution proceeds with the /-th instruction 
of the program concerned. If ##1 is itself the /-th instruction, then deadlock 
occurs. If / equals 0 or / is greater than the length of the program, then 
termination occurs. 

We define the meaning of PGLD programs by means of a function 
pgld2pga from the set of all PGLD programs to the set of all PGA programs. 
This function is defined by 


pgld2pga(ur;...; ue) = (Wir) ;...5 de(ur) 15)", 


where the auxiliary functions ~; from the set of all primitive instructions of 
PGLD to the set of all primitive instructions of PGA are defined as follows 
(Legh) 


HALL) = | ifl1=OVI>k, 


ua) te if wis not a jump instruction . 


Let P be a PGLD program. Then pgld2pga(P) represents the mean- 
ing of P as a PGA program. The intended behaviour of P is the be- 
haviour of pgld2pga(P). That is, the behaviour of P, written |Plpazp, 
is |pgld2pga(P)|. 

We use the phrase projection semantics to refer to the approach to 
semantics followed in this section. The meaning function pgld2pga is called 
a projection. 


5 Threads Interacting with Action Transforming 
Services 


A thread may perform certain actions for having itself affected by some 
service or affecting some service. When processing an action performed by 
a thread, an action transforming service affects that thread in one of the 
following ways: (i) by returning a reply value to the thread at completion 
of the processing of the action performed by the thread; (ii) by turning the 
processed action into another action. In this section, we introduce a mecha- 
nism that allows for services to affect threads in either way. The mechanism 
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is a generalization of the use mechanism introduced in [8].° We also intro- 
duce a complementary mechanism of that generalized use mechanism and an 
abstraction mechanism. The difference between the generalized use mech- 
anism and the complementary mechanism is a matter of perspective: the 
former is concerned with the effect of services on threads and therefore pro- 
duces threads, whereas the latter is concerned with the effect of threads on 
services and therefore produces services. The abstraction mechanism serves 
for concealment of the presence of internal actions, which arise primarily 
from the generalized use mechanism. 

We will use the generalized use mechanism and the complementary 
mechanism of that mechanism for interactions between threads that are 
behaviours exhibited by sequential programs on execution and services that 
process actions that relate to molecular dynamics. 

It is assumed that there is a fixed but arbitrary finite set of foct F and 
a fixed but arbitrary finite set of methods M. Each focus plays the role 
of a name of a service provided by the execution environment that can be 
requested to process a command. Each method plays the role of a command 
proper. For the set A of actions, we take the set {f.m| fe€ F,me My}. A 
thread performing an action f.m is considered to make a request to a service 
that is known to the thread under the name f to process command m. 

We introduce yet another sort: the sort S of services. S is considered 
to stand for the set of all services. We identify services with pairs (Hy, H2), 
where H,: M+ — {T,F,M,B} and Hyj:M* — Atau, satisfying the following 
conditions: 


Yne Me 
(Ja € M* + Hi(a~(m)) =M = Val € M* Hy(a! ~ (m)) ¢ {T,F}), 


Va € Mt,méeMe(Ai(a) =B > Hi(a~(m)) =B), 


Va € Mt «(Hi(a) 4M © Ho(a) = tau) © 


Let H be aservice, and let Hy and H2 be the unique functions such that 
H = (M1, H2). Then we write rf(H) and af(H) for H; and Hg, respectively. 


5In later papers, the original use mechanism is also called thread-service composition. 

°We write D* for the set of all finite sequences with elements from set D and Dt for 
the set of all non-empty finite sequences with elements from set D. We use the following 
notation for finite sequences: () for the empty sequence; (d) for the sequence having d as 
sole element; o ~ o’ for the concatenation of finite sequences o and o’. 


57 


Given a service H and a method m € M, the derived service of H after 
processing m, written ZH, is defined by rf (32H) (a) = rf(H)((m) ~a) 
and af (2;H)(a) = af (H)((m) ~a). 

A service H can be understood as follows: 


e if rf(H)((m)) =T, then the request to process m is accepted by the 
service, a positive reply is produced, m is turned into tau, and the 
service proceeds as aH ‘ 


e if rf(H)((m)) = F, then the request to process m is accepted by the 
service, a negative reply is produced, m is turned into tau, and the 
service proceeds as oH : 


e if rf(H)((m)) =M, then the request to process m is accepted by the 
service, no reply is produced, m is turned into af(H)((m)), and the 
service proceeds as =~ H; 


e if rf(H)((m)) = B, then the request to process m is rejected by the 
service. 


The three conditions imposed on services can be paraphrased as follows: 


e if it is possible that no reply is produced at completion of the process- 
ing of a command, then it is impossible that a positive or negative 
reply is produced at completion of the processing of that command; 


e after a request to process a command has been rejected, any request 
to process a command will be rejected; 


e areply is produced at completion of the processing of a command if 
and only if the command is turned into tau. 


We introduce the following additional constant and operators: 
e the divergent service constant D : S; 
e for each f € F, the binary use operator _ /f-:T x S — T; 
e for each f ¢ F, the binary apply operator _er_:TxS—S. 


We use infix notation for use and apply. 
D is the unique service H with the property that rf(H)(a) =B for all 
a € M*. The operators _ / — and _ ey _ are complementary. Intuitively, 
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Table 9: Axioms for use operators 


S/fH=S TSU1 
D/;H=D TSU2 
tauoa /; H =tauo(x/; H) TSU3 
(cdgmby)/sH = (e (EE Hf) if f#g TSU4 
(a di f.mb y) /¢ H =tauo (x /; 2 Bm Hf) if rf(H)((m)) =T TSU5 
(od fam y) /p H = tavoly /s 2H) if rf (H)((m)) =F TSUS 
(ed fom y) /p H = (w daf(H)((m)) © y) /p 2H if rf (H)((m)) =M_ TSUT 
@afimbey) /¢i =D if rf(H)((m)) =B TSU8 
(cdfmby)/,D=D TSU9 


Table 10: Axioms for apply operators 


Sep H =H TSA1 
De, H=D TSA2 
tauor)er, H=xe, H TSA3 
ridgmby)ey H=D if fAg TSA4 
cifme yep H=xey 3H if rf(H)((m)) =T TSAS5 
xifmey) ons 3H if rf(H)((m)) =F TSAG 
2 SfmB p) eH = Sof HMC) a) Belt Hf AY) = M TAT 
riaf.me y) e- if rf(H)((m)) =B TSA8 
cifmby)erD= D TSA9 
An>o (x) ef H=D) > ces H=D TSA10 


p /; H is the thread that results from processing all actions performed by 
thread p that are of the form f.m by the service H. When an action of 
the form f.m performed by thread p is processed by the service H, that 
action is turned into another action and, if this action is tau, postconditional 
composition is removed in favour of action prefixing on the basis of the reply 
value produced. Intuitively, pes H is the service that results from processing 
all basic actions performed by thread P that are of the form f.m by service 
AH. When an action of the form f.m performed by thread p is processed by 
service H, that service is changed into 3H : 

The axioms for the use and apply operators are given in Tables 9 and 10. 
In these tables, f and g stand for an arbitrary foci from F, m stands for an 
arbitrary method from M, and H is a variable of sort S. Axioms TSU3 and 


59 


Table 11: Axiom for abstraction 


Tiau(S) = S$ TT. 
fau(D) =D 112 
Trau(H <Itau y) = Trau(2) TT3 
Trau(@ Jal y) = tau(r) Jab tau(y) TT4 


TSU4 express that the action tau and actions of the form g.m, where f ¥ g, 
are not processed. Axioms TSU5—TSU7 express that a thread is affected by 
a service as described above when an action of the form f.m performed by 
the thread is processed by the service. Axiom TSU8 expresses that deadlock 
takes place when an action to be processed is not accepted. Axiom TSU9 
expresses that the divergent service does not accept any action. Axiom 
TSA8 expresses that a service is not affected by a thread when the action 
tau is performed by the thread and axiom TSA4 expresses that a service 
is turned into the divergent service when an action of the form g.m, where 
f #4 g, is performed by the thread. Axioms TSA5-TSA7 express that 
a service is affected by a thread as described above when an action of the 
form f.m performed by the thread is processed by the service. Axiom TSA8 
expresses that a service is turned into the divergent service when an action 
performed by the thread is not accepted. Axiom TSA9 expresses that the 
divergent service is not affected by a thread when an action of the form f.m 
is performed by the thread. 

We say that pe, H is a divergent thread application if t,(p) ef H = D 
for all n € N. Axiom TSA10 can be read as follows: if pe, H is a divergent 
thread application, then it equals D. 

The use operators introduced in [8] deals in essence with services H 
where af (H)(a) = tau for all a € M™. For these services, the use operators 
introduced here coincide with those use operators. 

Let T stand for either BTA, BTA+REC or BTA+REC+AIP. Then 
we will write T + TSI for T, taking the set {f.m | f € F,m € M} for A, 
extended with the divergent service constant, the use and apply operators, 
and the axioms from Tables 9 and 10, with the exception of axiom TSA10 
in the case where T does not stand for BTA+REC+AIP. 

The action tau is an internal action whose presence matters. To conceal 
its presence in the case where it does not matter after all, we also introduce 
the unary abstraction operator Tray :T — T. 

The axioms for the abstraction operator are given in Table 11. In this 
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table, a stands for an arbitrary basic action from A. 

Abstraction can for instance be appropriate in the case where tau arises 
from turning actions of an auxiliary nature into tau on use of a service. An 
example of this case will occur in Section 9. Unlike the use mechanisms 
introduced in [8] and in the current paper, the use mechanism introduced 
in [9] incorporates abstraction. 

Let T stand for either BTA, BTA+REC, BTA+REC+AIP, BTA+TSI, 
BTA+REC+TSI or BTA+REC+AIP+TSI. Then we will write T+ABSTR 
for T extended with the abstraction operator and the axioms from Table 11. 

The equation 7tay(tau”) = D is derivable from the axioms of BTA+ 
REC+AIP+ABSTR. 

To simplify matters, from now on the set {f.m | f € F,m € M} is 
taken as the set 21 of basic instructions when PGA or PGLD is concerned. 
Thereby no real restriction is imposed on the set 2: in the case where the 
cardinality of F equals 1, all basic instructions have the same focus and the 
set M of methods can be looked upon as the set 2 of basic instructions. 
We use the convention to omit foci from PGA programs in which all basic 
instructions have the same focus. 

Strictly speaking, the propositions and theorems presented in this pa- 
per are proved in the algebraic theory obtained by: (i) combining PGA with 
BTA+REC+AIP+TSI+ABSTR, resulting in a theory with three sorts: a 
sort P of programs, a sort T of threads, and a sort S of services; (ii) extend- 
ing the result by taking |_| for an additional operator from sort P to sort T 
and taking the semantic equations and rule defining thread extraction for 
additional axioms. 


6 State-Based Description of Services 


In this section, we introduce the state-based approach to describe families 
of services that will be used in Section 7. This approach is similar to the 
approach to describe state machines introduced in [9]. 

In this approach, a family of services is described by 


e aset of states S; 
e an effect function eff :M x S — S; 
e a yield function yld: M x S — {T,F,M,B}; 


e an action function act: M x S > Atau; 
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satisfying the following conditions: 


Ym € Me (As € S* yld(m,s) =M = Vs’ € S'« yld(m, s') ¢ {T, Fh) , 


dsEeSeVmeMe 
(yld(m,s) = BAYVs' € S'«(yld(m,s') =B => eff(m,s’) =s)), 


Ym € M,s € S«(yld(m,s) AM = act(m,s) = tau) . 


The set S contains the states in which the services may be, and the functions 
eff, yld and act give, for each method m and state s, the state, reply and 
action, respectively, that result from processing m in state s. 

We define, for each s € S, a cumulative effect function ceff,: M* —- S 
in terms of s and eff as follows: 


( 
ceff .(a o (m)) = eff (m, ceff ,(@)) ; 


We define, for each s € S, a service H, in terms of ceff,, yld and act as 
follows: 

rf(Hs)(a~ (m)) = yld(m, ceff ,(@)) , 

af (Hs)(a~ (m)) = act(m, ceff ,(@)) - 


H, is called the service with initial state s described by S, eff, yld and act. 
We say that {H, | s € S} is the family of services described by S, eff, yld 
and act. 

The conditions that are imposed on S, eff, yld and act imply that, for 
each s € S, H, is a service indeed. It is worth mentioning that 3H — 
Feg(m,s)> Tf(Hs)((m)) = yld(m, s), and af(Hs)((m)) = act(m, s). 


7 Services for Molecular Dynamics 


In this section, we describe a family of services which concerns molecular 
dynamics. The formal description given here elaborates on an informal 
description of molecular dynamics given in [1]. 

The states of molecular dynamics services resemble collections of 
molecules composed of atoms and the methods of molecular dynamics ser- 
vices transform the structure of molecules like in chemical reactions. An 
atom can have fields and each of those fields can contain an atom. An atom 
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together with the ones it has links to via fields can be viewed as a sub- 
molecule, and a submolecule that is not contained in a larger submolecule 
can be viewed as a molecule. Thus, the collection of molecules that make up 
a state can be viewed as a fluid. By means of methods, new atoms can be 
created, fields can be added to and removed from atoms, and the contents 
of fields of atoms can be examined and modified. A few methods use a 
spot to put an atom in or to get an atom from. By means of methods, the 
contents of spots can be compared and modified as well. Creating an atom 
is thought of as turning an element of a given set of proto-atoms into an 
atom. If there are no proto-atoms left, then atoms can no longer be created. 

It is assumed that a finite set Spot of spots such that F,M C Spot, 
a total order < on Spot, and a finite set Field of fields have been given. 
It is further assumed that a countable set PAtom of proto-atoms such that 
| ¢ PAtom and a bijection patom : [1,card(PAtom)] — PAtom have been 
given.’ Although the set of proto-atoms may be infinite, there exists at 
any time only a finite number of atoms. Each of those atoms has only a 
finite number of fields. Molecular dynamics services accept the following 
methods: 


e for each s € Spot, a create atom method s!}; 

e for each s,s’ € Spot, a set spot method s = s'; 

e for each s, € Spot, a clear spot method s = 0; 

e for each s,s’ € Spot, an equality test method s==s'; 

e for each s € Spot, an undefinedness test method s==0; 

e for each s € Spot and v € Field, a add field method s/v; 

e for each s € Spot and v € Field, a remove field method s\v; 

e for each s € Spot and v € Field, a has field method s|v; 

e for each s,s’ € Spot and v € Field, a set field method s.v = s’; 
e for each s,s’ € Spot and v € Field, a get field method s = s'.v; 


e for each s,s’ € Spot, a generate action method %s.%s'. 


"We use the notation [n,m], where n,m €N, for the set {i EN|n <i < my}. 
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We write Mgen for the set of all generate action methods and Mmg for 
the set of all methods of molecular dynamics services. It is assumed that 


Mmd Cc M. 


The states of molecular dynamics services comprise the contents of 


all spots, the fields of the existing atoms, and the contents of those fields. 
The methods accepted by molecular dynamics services can be explained as 
follows: 


s!: if an atom can be created, then the contents of spot s becomes 
a newly created atom and the reply is T; otherwise, nothing changes 
and the reply is F; 


s = s': the contents of spot s becomes the same as the contents of 
spot s’ and the reply is T; 


s =0: the contents of spot s becomes undefined and the reply is T; 


s==s': if the contents of spot s equals the contents of spot s’, then 
nothing changes and the reply is T; otherwise, nothing changes and 
the reply is F; 


s==0: if the contents of spot s is undefined, then nothing changes 
and the reply is T; otherwise, nothing changes and the reply is F; 


s/v: if the contents of spot s is an atom and v is not yet a field of that 
atom, then v is added (with undefined contents) to the fields of that 
atom and the reply is T; otherwise, nothing changes and the reply is 
F; 


s\v: if the contents of spot s is an atom and v is a field of that atom, 
then v is removed from the fields of that atom and the reply is T; 
otherwise, nothing changes and the reply is F; 


s|v: if the contents of spot s is an atom and v is a field of that atom, 
then nothing changes and the reply is T; otherwise, nothing changes 
and the reply is F; 


s.v = s‘: if the contents of spot s is an atom and v is a field of that 
atom, then the contents of that field becomes the same as the contents 
of spot s’ and the reply is T; otherwise, nothing changes and the reply 
is F; 
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e s=s'.v: if the contents of spot s’ is an atom and v is a field of that 
atom, then the contents of spot s becomes the same as the contents 
of that field and the reply is T; otherwise, nothing changes and the 
reply is F; 


e %s.%s': if the contents of spots s and s’ are atoms and there exist 
f' € F and m’ € M such that the contents of spot f’ equals the 
contents of spot s and the contents of spot m’ equals the contents of 
spot s’, then nothing changes, there is no reply, and %s.%s’ is turned 
into f.m where f and m are the least f’ € F and m’ € M with respect 
to < that satisfy the conditions formulated above; otherwise, %s.%s' 
is rejected. 

In the explanation given above, wherever we say that the reply is T or the 
reply is F, this is meant to imply that the method concerned is turned into 
tau. Moreover, wherever we say that the contents of a spot or field becomes 
the same as the contents of another spot or field, this is meant to imply that 

the former contents becomes undefined if the latter contents is undefined. 

Let 


SS  =Spot > (PAtomU{L}), 


AS = |) (4> LU (F= (PAtomu {1}))), 


A€P xn (PAtom) F EP ¥n (Field) 
MDS = {(0,a) € SS x AS | rng(o) C dom(a) U{L} A 
Va € dom(a) «rng(a(a)) C dom(a) U{L}}, 


s € MDS. 


Then we write MD, for the service with initial state s described by S = 
MDS U {tT}, where | ¢ MDS, and the functions eff, yld and act defined in 
Tables 12, 13 and 14.° In these tables, the functions new : Pan(PAtom) > 
(PAtomU{1}), gencnd: SS x Spot x Spot > {T,F}, and genact: SS x Spot x 
Spot — A are used. These functions are defined as follows: 


new(A) = patom(m + 1) if m < card(PAtom) , 
new(A) = 1 if m > card(PAtom) , 


’We use the following notation for functions: [] for the empty function; [d +> r] for 
the function f with dom(f) = {d} such that f(d) = r; f @g for the function h with 
dom(h) = dom(f) Udom(g) such that for all d€ dom(h), h(d) = f(d) if d ¢ dom(g) and 
h(d) = g(d) otherwise; and f <D for the function g with dom(g) = dom(f) \ D such that 
for all d € dom(g), g(d) = f(d). 
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Table 12: Effect function for molecular dynamics services 


eff (s!,(0,a)) = 

(o @ [s & new(dom(a))], a © [new(dom(a)) + []]) if new(dom(a)) A L 
eff (s!, (7, a)) = (7, a) if new(dom(a)) = L 
eff (s = s',(0,a)) = (o 6 [s & o(s’)], a) 
eff (s = 0, (0,a)) = (0 @ [s & L], a) 
eff (s==s', (0, a)) = (0, a) 
eff (s==0, (0, a)) = (g,a) 
eff (s/v, (a, a)) = 

(2,06 [o(s) + a(o(s)) @ [v > LI) if o(s) £ LAv ¢ dom(a(o(s))) 
eff (s/v, (a, a)) = (0, a) if o(s) = LV v € dom(a(a(s))) 
eff (s\v, (a, a)) = (0,a @ [a(s) & a(o(s)) = {v}]) if o(s) # L Av € dom(a(a(s))) 
eff (s\v, (a, @)) = (0, a) if o(s) = LV v ¢ dom(a(o(s))) 
eff (s|v, (0, @)) = (a, a) 
eff (s.v = s',(o,a)) = 

(c, 0.6 [o(s) + a(o(s)) @ [v o(s")]) if o(s) # LAv € dom(a(o(s))) 
eff (s.v = 8’, (0, a)) = (a, a) if o(s) = LV v ¢ dom(a(o(s))) 
eff (s = s'.v, (0, a)) = (6 6 [s & a(o(s’))(v)], @) if o(s') A L Av € dom(a(a(s’))) 
eff (s = s’.v, (0, a)) = (0, a) if o(s') = LV v €¢ dom(a(a(s’))) 
eff (%s.%s' , (o, w)) = (0, a) if gencnd(o, s,s’) = T 
eff (%s.%s', (o,a)) = T if gencnd(o, s,s’) =F 
eff (m, (a, a)) = TF if m € Mma 
eff(m, 1) = 


where m = max{n | patom(n) € A}; 


gencnd(a, s,s’) = T iff 
a(s) # LAo(s') A LAA € Feo(f) =a(s) Adm € Mea(m) = a(s’) ; 


genact(o, s,s’) = f.m, 


where f = min{f’ € F | o(f’) = o(s)} and m = min{m’ € M | o(m’) = 
o(s)}. We write MD for the family of services described by S, eff, yld, and 
act. We write MD, for the service from this family of which the initial 
state is the unique (¢,a) € S such that dom(a) = []. 

Let (o,a) € S, let s € Spot, let a € dom(a), and let v € dom(a(a)). 
Then o(s) is the contents of spot s if o(s) 4 L, v is a field of atom a, and 
a(a)(v) is the contents of field v of atom a if a(a)(v) # L. The contents 
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Table 13: Yield function for molecular dynamics services 


yld(s!,(o,a)) =T if new(dom(a)) A 
yld(s!,(o,a)) =F if new(dom(a)) = 

yld(s = s’,(o,a)) =T 

yld(s = 0, (o,a)) =T 

yld(s==s',(o,a))=T if o(s) = 0(s’) 
yld(s==s',(o,a))=F_ if a(s) £a(s’) 
yld(s==0,(0,a))=T if o(s)= +L 
yld(s==0,(0,a))=F if o(s) 4 4 

yld(s/v,(o,a)) =T if o(s) A LAv ¢ dom(a(o(s))) 
yld(s/v,(o,a)) =F if o(s) = L Vv € dom(a(a(s))) 
yld(s\v, (o,a)) =T if o(s) # L Av € dom(a(o(s))) 
yld(s\v, (o,a)) =F if o(s) = LV v ¢ dom(a(o(s))) 
yld(s|v, (o,a@)) =T if o(s) # L Av € dom(a(a(s))) 
yld(s|v, (o,a)) =F if o(s) = LV v ¢ dom(a(o(s))) 
yld(s.v = s’,(0,a)) =T if o(s) A LAv € dom(a(a(s))) 
yld(s.v = s’,(0,a)) =F if o(s) = 1 Vv ¢ dom(a(o(s))) 
yld(s = s’.v,(o,a)) =T if o(s’) 4 L Av € dom(a(a(s’))) 
yld(s = s’.v,(o,a)) =F if o(s’) = LV v ¢ dom(a(a(s’))) 
yld(%s.%s',(o,a)) =M if gencnd(o, s,s’) = T 
yld(%s.%s',(o,a)) =B_ if gencnd(o, s,s’) = F 

yld(m, (o,a)) = B ifm ¢ Mma 

yld(m,t) =B 


Table 14: Action function for molecular dynamics services 


act(m, (a, a)) = tau if m € Meen 
act(%s.%s', (a, v)) = genact(o, s,s’) if gencnd(c, s,s’) =T 
act(%s.%s' , (o, w)) = tau if gencend(o, s,s’) = F 
act(m, 7) = tau 


of spot s is undefined if o(s) = L, and the contents of field v of atom a 
is undefined if a(a)(v) = L. Notice that dom(qa) is taken for the set of 
all existing atoms. Therefore, the contents of each spot, i.e. each element 
of rng(o), must be in dom(a) if the contents is defined. Moreover, for 
each existing atom a, the contents of each of its fields, i.e. each element of 
rng(a(a)), must be in dom(q) if the contents is defined. Molecular dynamics 
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services get into state T when refusing a request to process a command. The 
molecular dynamics service with initial state | is the same as the divergent 
service D. The function new turns proto-atoms into atoms. After all proto- 
atoms have been turned into atoms, new yields L. This can only happen 
if the number of proto-atoms is finite. The initial state of MDjnj is the 
unique state in which no proto-atoms have been turned into atoms yet. 


8 Representing PGLD Programs by Molecules 


In this section, we associate each PGLD program with a PGA program for 
constructing a representation of the PGLD program by a molecule. 

Let uz ;...3; uz be a PGLD program in which the foci f1,..., f, and 
the methods m,,...,™,y occur. The idea is that: 


e an atom is created for each of the foci fi,...,f, and the methods 
My ,..., Mn’, using the focus or method itself as the spot into which 
the atom concerned is brought on its creation, and an atom is created 
for each of the instructions uj1,..., ug; 


e each atom that corresponds to an instruction is linked via fields to 
other atoms as follows: 


— if the corresponding instruction is of the form f.m, +f.m or 
—f.m, then the atom is linked to the atoms that correspond to 
the focus and method concerned and the atoms that correspond 
to the instructions with which execution must proceed in the 
cases of a positive and a negative reply; 


— if the corresponding instruction is of the form ##/ and 1 < 
1 <k, then the atom is linked to the atom that corresponds to 
the instruction with which execution must proceed, i.e. the /-th 
instruction; 


e to prevent that an exception must be made of the instruction wu, in 
the case where it is of the form f.m, +f.m or —f.m, two additional 
atoms are created that are not linked to other atoms; 


e the atom that corresponds to the first instruction is brought into spot 
8. 


Notice that, if an atom corresponds to an instruction of the form ##/ and 
not 1<1<k, then the atom is not linked to another atom. 
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We make the assumptions that F and M are disjoint and that PAtom is 
infinite. Under these assumptions, atom creation always leads to a positive 
reply and no test instructions are needed in the programs for constructing 
representations of PGLD programs. The first assumption is only made 
because it permits this simplification. The second assumption is primarily 
made because there will be PGLD programs for which there are no PGA 
programs for constructing their representation if PAtom is finite. 

We define the PGA programs for constructing representations of PGLD 
programs by means of a function pgld2md from the set of all PGLD programs 
to the set of all PGA programs. This function is defined by 


Pebdsmd (ays ohio tp SH Fe eat ee ie see eT BTS, aes 
pi(ur) ;--+5 Pe(ux) 5 
8p41/ stop ; 8z42/ stop ; 
s=51;!, 


where the auxiliary functions p; from the set of all primitive instructions of 
PGLD to the set of all PGA programs are defined as follows (1 < j < k): 


pi(f.m) = s8;/focus ;s;/meth ; s;/pos ; s;/neg ; 
$;.focus = f ;s;.meth = ™ ; 8;.pos = 8341 5 8;.neg = $541, 
pj(+f.m) = s;/focus ; s;/meth ; s;/pos ; s;/neg ; 
8;.focus = f ; s;.meth = ™ ; 8;.pos = 8541 5 8;.neg = 8542, 
pj(—f.m) = s;/focus ; s;/meth ; s;/pos ; s;/neg ; 
8;.focus = f ;s;.meth = ™ ; 8;.pos = 8j12 5 8;.neg = 8741, 


pi(##HL) = s;/jmp;sj.jmp=s ifl<Il<k, 


pi(##1) = 5;/stop if “(1 <1<k) 
and 
e fi,...,fn © F are the different foci that occur in uy ;...; Ug; 
© m1,..-,Mn € M are the different methods that occur in ui ;...3; ug; 


@ 8,51,...,Sk42 € Spot \(FUM). 


In this definition, the omitted focus is considered to be md. 

The properties stated in the following proposition show that it is easy 
to retrieve the PGLD program P from its representation constructed by 
pgld2md(P). 
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Proposition 1 Let P = u,;...; ug be a PGLD program, and let (a, a) be 
such that MD (¢.4) = |pgld2md(P)|¢mgMDinit- Then yld(s==s1,(o,a)) =T 
and for all j,l€ [1,k], fe F andme mM: 


eu; = f.m iff 
yld(u = s;.focus,(o,a)) =T, yld(u==f, eff(u = s;.focus,(o,a))) =T 
yld(v = s;.meth, (o,a)) =T, yld(v==m, eff (u = s;.meth, (o,a))) =T, 
yld(s = s;.pos,(o,a))=T, yld(s==s,41, eff (s = s;.pos, (0, a)) 
yld(s = s;.neg,(o0,a))=T, yld(s==s,41, eff (s = s;.neg, (o,a))) = 


eu; =+f.m iff 
yld(w= s;.focus,(a,a)) = T, yld(u==f, eff (u = s;.focus,(o,a))) =T 
yld(v = s;.meth,(o,a)) =T, yld(v==m, eff (u = s;.meth, (o,a))) =T 
yld(s = s;.pos,(o,a)) = 

( 


T,  yld(s==s;+1, eff (s = 8;.pos,(o,a))) =T, 
yld(s = s;.neg,(o0,a))=T, yld(s==s;+2, eff (s = s;.neg, (o,a))) =T; 
e Uj = —f.m iff 
yld(u = s;.focus,(o,a)) =T, yld(u==f, eff(u = s;.focus,(o,a))) =T, 
yld(v = s;.meth, (o,a)) =T, yld(v==m, eff (u = s;.meth, (o,a))) =T, 
yld(s = s;.pos,(o,a))=T, yld(s==s;+2, eff(s = s;.pos,(o,a))) =T, 
yld(s = 8;.neg,(o,a))=T, yld(s==s;41, ef (8 = s8;.neg, (0, a))) = T; 


eu; = ##l iff 
yld(s = s;.jmp,(o,a)) =T, yld(s==s), eff (s = s;.jmp, (o,))) = T; 


e uj =##l' for some I’ ¢ [1,k] iff yld(s;|stop, (o,a)) = T. 


Proof: From the assumption that F and M are disjoint, the assumption 
that PAtom is infinite, and the definition of pg1d2md, it follows easily that: 


e all atom creations are successful; 


e the content of each spot used in an atom creation does not change 
after the atom creation concerned; 


e p;(u;) modifies the atom that is the contents of spot s; only (1 <j < 
k), 8p41/stop modifies the atom that is the contents of spot 5,41 only, 
and s,19/stop modifies the atom that is the contents of spot 449 only. 


From this, the properties stated in the proposition follow straightforwardly 
by case distinction as in the definition of p;. 


70 


9 Interpreting PGLD Programs Represented by 
Molecules 


In this section, we introduce a PGA program for interpreting PGLD pro- 
grams represented by molecules. 
The idea is that: 


e if the atom to be interpreted has the field stop, then the interpretation 
terminates; 


e if the atom to be interpreted has the field jmp, then the interpretation 
proceeds with the atom to which it is linked via this field; 


e otherwise, first the basic instruction represented by the atoms to which 
it is linked via the fields focus and meth is executed and then the 
interpretation proceeds with the atom to which it is linked via the 
field pos or the field neg, depending upon the reply being positive or 
negative. 


The following is the PGA program for interpreting PGLD programs 
represented by molecules: 


(+8|stop ;!; 
+s|jmp ; #9 ; 
u = 8.focus ;v = s.meth;+%u.%v ; #3; 5 = s.neg; #4; 5 = s.pos ; #2; 


8 = s.jmp)* , 


where u,v € Spot \(F UM). Again, the omitted focus is considered to be 
md. Below, we write I for this PGA program. 

The program I interprets PGLD programs correctly in the sense that, 
for all PGLD programs P, the intended behaviour of P under execution 
coincides with the behaviour of the interpreter under execution on inter- 
action with the molecular dynamics service that holds the representation 
of P constructed by the program pgld2md(P) when abstracted from tau. 
This is stated rigorously by Theorem 1 below. The theorem is preceded by 
Proposition 2 which is used in the proof of the theorem. The proposition 
concerns the local use of the spots u and v in I. 

For convenience, we introduce some abbreviations. Let P = u,;...; 
ux be a PGLD program and let i € [1,k]. Then we write MD,.p,(p) for 
|pgld2md(P)| @ma MDinit and MDyepr(;,p) for |s = $i; !| (ma MD repr(p)- 
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Proposition 2 Let P= ,;...; uz be a PGLD program, let i € [1,k], and 
let MD" opx(i,P) be such that MDyepr(i,p) = |U=0;0 =0;!| ma MD! opx(i,P): 


Then Trau(|Z| /md MD yepr(i,P)) = Trau(|Z| /ma MD epr(i,P)): 


Proof: By AIP, it is sufficient to prove that for all n € N, for all 7 € [1, A], 


Tn (Ttau(|Z| /md MD repr(i,P))) = Tnr(Ttau(|Z| /ma MD epr(i,P))): This property 
is straightforward to prove by induction on n. 


Theorem 1 For all PGLD programs P: 
|Plpcrp = Trau(|L| /md MD yepr(P)) : 


Proof: In the proof, we make use of an auxiliary function |_,_| from the 
set of all natural numbers and the set of all PGLD programs to the set of 
all closed terms of sort T. It gives, for each natural number 7 and PGLD 
program u,;...; Ug, a closed term of sort T that denotes the behaviour of 
U1;...}Uk when executed from position 7 if 1 <i < k and S otherwise. This 
function is defined as follows: 


[eee So Gti | eG) eens ees els Gi Gi) sees i) else | 
ifl<i<k, 
iit UES if 7»l<i<k 


(where ~; is as in the definition of pgld2pga). It follows easily from the 
definition of |_,_| and the axioms of PGA that if 1 <i<k: 


a, thy genes oe| So t+ Ty urs ss peg ifuj=a, 

Ji,ur;..-3Up| = |e +1,ur;...; ug] dab |t+2,u,;...; uz) if u=t+a, 

Ji,Ur;..-3 UR| = |@+2,u,;...5 ug] Jab t+ 1,u.;...5; ug) if u=—a, 

Petia Sos. 2s tip ae ns ae if u = ##1. 
Let P = uj, ;...; uz be a PGLD program, let 


T = {|i, P| [ie [1 A]}, 
T'= {Ttau (|Z| /md MD yepr(i,P)) tE Ls kl} ) 


and let @:T — T’ be the bijection defined by 
B( li, P|) = Tau (|Z| /md MD yepr(i,P)) : 
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For each closed term p’ of sort T, write 6*(p’) for p’ with, for all p € T, all 
occurrences of p in p’ replaced by 3(p). Then, using the equations concerning 
the auxiliary function |_,_| given above and Propositions 1 and 2, it is 
straightforward to prove that there exists a set E consisting of one derivable 
equation p = p’ for each p € T such that, for all equations p = p’ in E: 


e the equation 3(p) = B*(p’) is derivable; 


e p' €T only if p’ can be rewritten to a p’ ¢ T using the equations in 
FE from left to right. 


Because |Pleaxp = 1, P| and Trau (|Z| /md MD yepr(P)) =; Trau(|Z| /md 
MD yepr(1,P))s this means that |P|pcrp and Ttau(|L| /ma MD yepr(P)) are so- 
lutions of the same guarded recursive specification. Because guarded re- 
cursive specifications have unique solutions, it follows immediately that 
|P|peup = Trau (|| /md MD yepr(P))- 


In the program pgld2md(u; ;...; uz), we could have replaced s = s; ;! 
by s = 8135, =03...3 5442 = 0;! because the program I does not use 
the spots 51,..., 542. However, we cannot conceive a proof of Theorem 1 


in which these auxiliary spots are not found. The main problem is that 
we cannot find a way of formulating the gist of Proposition 1, which looks 
to be material to a proof of the theorem, without referring to the spots 
S1,+++,S5k+2- 


10 PGLD with Indirect Jumps 


In this section, we introduce a variant of PGLD with indirect jump instruc- 
tions. This variant is called PGLD;;. However, preceding the introduction 
of PGLD,;, we give a state-based description of the very simple family of 
services that constitute a register file of which the registers can contain nat- 
ural numbers up to some bound. This register file will be used later on to 
describe the behaviour of PGLDj; programs. 

It is assumed that a fixed but arbitrary number J has been given, 
which is considered the number of registers available. It is also assumed 
that a fixed but arbitrary number N has been given, which is considered 
the greatest natural number that can be contained in a register. 

The register file services accept the following methods: 


e for each i € [0,7] and n € [0, NJ, a register set method set:i:n; 
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e for each i € [0,7] and n € [0, N], a register test method eq:i:n. 


We write Mregs for the set {set:i:n,eq:i:n | 7 € [0,I] An € [0,N]}. It is 
assumed that Mregs C M. 

The methods accepted by register file services can be explained as fol- 
lows: 


e set:i:n: the contents of register 1 becomes n and the reply is T; 


e eq:i:n: if the contents of register 7 equals n, then nothing changes and 
the reply is T; otherwise nothing changes and the reply is F. 


Let s:{1,Z] — [0, NJ]. Then we write Regs, for the service with initial 
state s described by S = ({1,J] — [0,N]) U{t}, where f ¢ [1,1] — [0, N], 
and the functions eff, yld and act defined as follows (i € [0,Z], n € [0, N], 
p: (1, I] — [0, NJ): 


eff (set:i:n, p) =p@lirn], 


( 
eff (eq:i:n, p) =p, 
eff (m, p) = T ifm ¢ Mregs 5 
eff(m,T) =T, 
yld(set:i:n, p) =T , 
yld(eq:i:n, p) = T if p(i) =n, 
yld(eq:i:n, p) = F if la) es 
yld(m, p) = B ifm ¢ Mregs ; 
yld(m, 1) =B, 
act(m, p) = tau , 
act(m, t) = tau 


We write Regsini, for Regsi1_.o}6...e [140]: 

We pass on to PGLD with indirect jump instructions. In PGLD,;, 
it is assumed that there is a fixed but arbitrary finite set of foci F with 
regs © F and a fixed but arbitrary finite set of methods M. Moreover, 
we adopt the assumptions made about register file services above. The set 
{f.m| f © F \ {regs}, m € M} is taken as the set 2% of basic instructions. 

PGLD;j; has the following primitive instructions: 


e for each a € 2, a plain basic instruction a; 
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e for each a € A, a positive test instruction +a; 


e for each a € A, a negative test instruction —a; 


for each 1 € N, a direct absolute jump instruction ##1; 
e for each i € [1,7] and n € [1, N], a register set instruction set:i:n; 
e for each i € [1, I], an indirect absolute jump instruction i###1. 


PGLD;j; programs have the form uj ;...; uz, where uji,...,uz are primitive 
instructions of PGLD). 

The plain basic instructions, the positive test instructions, the negative 
test instructions, and the direct absolute jump instructions are as in PGLD. 
The effect of a register set instruction set:i:n is that the content of register 
i becomes n. The effect of an indirect absolute jump instruction i##i is 
that execution proceeds with the /-th instruction of the program concerned, 
where | is the content of register i. If i###7 is itself the /-th instruction, 
then deadlock occurs. If | equals 0 or / is greater than the length of the 
program, termination occurs. We stipulate that the content of each register 
is initially 0. 

Like before, we define the meaning of PGLD,; programs by means of a 
function pgldij2pgld from the set of all PGLDj; programs to the set of all 
PGLD programs. This function is defined by 


pgldij2pgld(uy;...; upg) = 
(ui) ;...; blue) ; FO ; FH ; 
+regs.eq:1:1 ; #341;...; +regs.eq:1:n ; #4tn ; ##0 ; 
+regs.eq:[:1; ##1;...; +regs.eq:I:n ; ##n;##0, 


where n = min(k, NV) and the auxiliary function w from the set of all primi- 
tive instructions of PGLDj; to the set of all primitive instructions of PGLD 
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is defined as follows: 


(a) =a, 

w(ta) =+a, 

p(-a)  =-a, 

V(H#) = ##I flsk, 
OH!) = #0 ifl>k, 
w(set:i:n) = regs.set:i:n , 
vit) = ##L , 


and for each i € [1, J]: 
Ll =k+3+(2-min(k,N)+1)-(@-1). 


The idea is that each indirect absolute jump can be replaced by a direct 
absolute jump to the beginning of the instruction sequence 


+regs.eq:i:1 ; ##1;...; +regs.eq:i:n ; ##n;##0, 


where 7 is the register concerned and n = min(k, N). The execution of this 
instruction sequence leads to the intended jump after the content of the 
register concerned has been found by a linear search. To enforce termination 
of the program after execution of its last instruction if the last instruction 
is a plain basic instruction, a positive test instruction or a negative test 
instruction, ##0 ; #70 is appended to W(u1);...; W(ug). Because the 
length of the translated program is greater than k, care is taken that there 
are no direct absolute jumps to instructions with a position greater than k. 
Obviously, the linear search for the content of a register can be replaced by 
a binary search. 

Let P be a PGLD; program. Then pgldij2pgld(P) represents the 
meaning of P as a PGLD program. The intended behaviour of P is 
the behaviour of pgldij2pgld(P) on interaction with a register file when 
abstracted from tau. That is, the behaviour of P, written |Plraup,;+ is 
Trau([PEldij2pgld(P)|pcrp /regs Regsinit)- 

A slightly different variant of PGLD with indirect jump instructions is 
introduced in [4] under the same name. 
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11 The Interpretation of PGLD;; Programs 


In this section, we associate each PGLDj; program with a PGA program 
for constructing a representation of the PGLD;; program by a molecule and 
introduce a PGA program for interpreting PGLD;; programs represented 
by molecules. This amounts to enhancing the PGA programs given in Sec- 
tions 8 and 9. PGLD;; programs without indirect absolute jump instructions 
are represented and interpreted exactly as before. 

The idea of the enhancements is that: 


e an atom is created for each of the registers; 


e each atom that corresponds to a register is handled as if it concerns 
a direct jump instruction to the instruction with the content of the 
register as position; 


e each atom that corresponds to an indirect jump instruction is han- 
dled as if it concerns a direct jump instruction to the direct jump 
instruction that takes the place of the register concerned; 


e each register set instruction is handled as if it concerns an instruction 
for overwriting the direct jump instruction that takes the place of the 
register concerned. 


This means that the program for interpreting PGLDj; programs represented 
by molecules will change the molecular representation of the PGLDj; pro- 
gram being interpreted instead of the state of a register file whenever it 
comes across a register set instruction. 

We define the PGA programs for constructing representations of 
PGLD;; programs by means of a function pgldij2md from the set of all 
PGLDj; programs to the set of all PGA programs. This function is defined 
by 

pgldij2md(uz;...; uz) = 
Fite ees tol ess ig 2 Syn Spal ol sn es 
pi(ur) 3 ..-3 P(e) ; 
8p41/ Stop ; Sp42/stop ; 
s\/jmp 3 s.jmp = Sk425--.5 87/Jmp 5 $7.JMp = SK42 5 


s=581;}, 
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where the auxiliary functions p; from the set of all primitive instructions of 
PGLDj; to the set of all PGA programs are defined as follows (1 < j < k): 


pj(f.m) = 8;/focus ;s;/meth ; s;/pos ; s;/neg ; 
8;.focus = f ; 8;.meth = ™ ; 8;.pos = 8341 3 8j.neg = 854 
pj(+f.m) = s;/focus ; s;/meth ; s;/pos ; s;/neg ; 


s;.focus = f ; 8;.meth =m; 8}.pos = 8j41 3 8j.neg = 8j4 


oe 
pj(—f.m) = s;/focus ; s;/meth ; s;/pos ; s;/neg ; 

8;.focus = f ; 8;.meth = ™ ; 8;.pos = 8j49 ; 8j.neg = $741 , 
pil) = 8;/Jmp ; 8;.Jmp = 8 ss Se an 
pi(##l) = 5;/stop if SR) 


p;(set:i:l) = s;/reg ; s;/cont ; s;/nat ; 

Sp.feg = 3, 5,00 = 35s net = 87 HPL 
p;(set:i:l) = s;/reg ; s;/cont ; s;/nat ; 

$).7eg- = 8, 38;.ont = S195 sent Sepa if a(S < ky, 
pi (i##i) = sj /jmp ; s;.jmp = 8; 


and 
e fi,...,fn © F are the different foci that occur in uz ;...; ug; 
© ™1,...,Mn € M are the different methods that occur in uj ;...3 ug; 


© S,S1,..-,Sk42,$4,---,8) € Spot \(FUM). 


The following is the PGA program for interpreting PGLD;; programs 
represented by molecules: 


(+s|stop ;!; 

+s|jmp ; #16 ; 

+s|reg ; #9 ; 

u = s.focus;v = s.meth;+%u.%vu ; #338 = s.neg;#9;s =s.pos; #7; 


uUu=s.reg;v = s.cont;ujmp =v;s = s.nat; #2; 


s = s.jmp)” , 


where u,v € Spot \(F UM). Below, we write I’ for this PGA program. 
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Theorem 2 below states rigorously that program J’ interprets PGLD;; 
programs correctly. In that theorem, other than in Theorem 1, we write 
MD yepr(p) for |pgldij2md(P)| ¢ma MDinit- 


Theorem 2 For all PGLDi programs P: 


|Plraup, = Trau(|L"| /md MD yepr(P)) : 


Proof: The proof follows the same line as the proof of Theorem 1. This 
is possible because on interpretation any change of the state of the register 
file is reflected by a corresponding change of its molecular representation. 


12 Conclusions 


In this paper, we have considered the programming of an interpreter for a 
program notation that is close to existing assembly languages using PGA 
with the primitives of molecular dynamics as basic instructions. We have 
given PGA programs for constructing representations of the programs to be 
interpreted by molecules and a PGA program for interpreting those repre- 
sentations and we have shown that the latter PGA program does the inter- 
pretation correctly. We have experienced that, although primarily meant 
for explaining programming language features relating to the use of dynamic 
data structures, the collection of primitives of molecular dynamics in itself 
is suited to the programming wants concerned. 

We observe that: (i) the program notation for which the presented 
interpreter has been designed belongs to the simplest program notations 
devised ever, (ii) it is hard to imagine that the programs to be interpreted 
can be represented by molecules in a way that is simpler than the way chosen 
for the presented interpreter, (iii) it is hard to conceive of an interpreter that 
is simpler than the presented interpreter. This means not at all that the 
design of the interpreter was simple. On the contrary, the design turned out 
to be disappointingly difficult. It happened that, owing to the quest for a 
simple interpreter, it was inescapable that the design was to a great extent 
a trial-and-error matter. 

Dynamic data structures modelled using molecular dynamics can 
straightforwardly be implemented in programming languages ranging from 
PASCAL [15] to C# [11] through pointers or references, provided that fields 
are not added or removed dynamically. Using molecular dynamics, we need 
not be aware of the existence of the pointers used for linking data. The 
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name molecular dynamics refers to the molecule metaphor used in the in- 
troduction. By that, there is no clue in the name itself to what it stands for. 
To remedy this defect, we suggest data linkage dynamics as an alternative 
name. 
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